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Hallen type integral equations are formulated for an arbitrarily 
bent monopole above an infinite perfect ground plane and* by numerical 
solution techniques* current distributions are obtained. The coupled 
integral equations are founded upon the familiar electric scalar poten- 
tial and magnetic vector potential. Numerical results are presented 
for several cases including different bent angles and are compared where 
possible with available data. 
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I. INTRODUCTION 



A linear antenna is one in which a "long dimension" can be discerned; 
for example, a long circular cylinder the length of which is much greater 
than the diameter is looked upon as a linear antenna. Initially, the em- 
phasis will be limited to the so-called cylindrical antenna driven at its 
mid-point by a voltage source V applied to a vanishingly short gap at the 
antenna center. The antenna will be considered to reside in a lossless 

> 

homogeneous and isotropic medium of infinite extent. 

A thin wire antenna is commonly defined as a linear antenna which has 
a radius much less than a wavelength. There are several approximations 
which may be employed in the analysis of thin wire antennas. 

(1) The current on the cylindrical surface is assumed to be 
entirely in the axial direction. 

(2) The current distribution is assumed to be constant around 
the periphery of the wire at a given cross-section and hence can be 
approximated as a filament of current on the axis of the cylinder. 

(3) The boundary conditions are applied only to those components 
of the total electric field which are tangent to the wire surface and 
parallel to the wire axis. 

Many investigators have considered the theoretical problem of deter- 
mining the current induced on a thin wire scattering antenna. [Ref. 1,2,3]. 
These investigations have been extended to include the treatment of ar- 
bitrary configurations of thin wire scattering antennas intersecting at 
a single junction. [Ref. 4]. Three particular configurations for which 
numerical results have been obtained are a set of nonintersecting thin 
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straight wires, [Ref. 5], a set of perpendicular intersecting thin straight 
wires, [Ref. 6], and a pair of skewed crossed wires. [Ref. 7]. 

One general approach to finding the current on thin wire antennas is 
to obtain a set of integral equations which involve the unknown current 
distributions. Using the equation of continuity, one may determine ex- 
pressions for the vector and scalar potentials in terms of the current 
distribution. It is in obtaining these integral equations that a choice 
can be made. 

Following the Pockington equation formulation, the expressions for 
the vector and scalar potentials are used along with Maxwell’s equations 
to derive a system of integro-dif ferential equations for the unknown 
current distribution. An auxiliary scalar function is defined and intro- 
duced into the system of integro-dif ferential equations. The solution of 
these now ordinary differential equations subject to appropriate boundary 
conditions gives the desired set of * integral equations for the unknown 
current distribution. 

If the Halldn integral equation approach is used, one independently 
establishes two expressions for a given component of the vector potential 
in the direction of a wire element and these expressions are evaluated on 
the wire surface. One expression is obtained from the solution of an 
inhomogeneous differential equation satisfied by the component of vector 
potential. 

Another expression for vector potential is obtained directly from the 
familiar vector potential itself. These two expressions must be equal 
and, hence, are equated. Next, the boundary condition that the tangential 
component of electric field on the surface of the wire element must be 
zero is applied. Now the arbitrary constants of integration which arise 
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in the solution of the inhomogeneous differential equation are evaluated. 
The constants of integration must have assumed values which are in com- 
pliance with all other required conditions. [Ref. 8], 

There have been numerous methods proposed for solving the integral 
equations, however obtained. These methods have fallen into a few broad 
catagories : 

(1) Iteration [Ref. 9^ 10]. 

(2) Expansion of the current in a Fourier series [Ref. 11]. 

(3) The variational method [Ref. 12,13,14]. 

(4) The Wiener-Hopf techniques for the case of the long 
antenna [Ref. 15]. 

(5) Numerical methods [Ref 1, 16] . 

One effective numerical method of solving the integral equation is 
the so-called direct integration technique (Method of Moments). [Ref. 16, 
17]. Basically the technique utilizes an expansion of the unknown func- 
tion in some convenient manner. This expression may be a finite sum of 
piecewise continuous functions which are non-zero over some small inter- 
val of definition and vanish elsewhere. A substitution of this expres- 
sion for the unknown function into the system of integral equations 
reduces them to a system of linear algebraic equations for the unknown 
coefficients. These linear algebraic equations may be solved by conven- 
tional techniques . 

The purpose of this thesis is to formulate the coupled integral 
equations for an arbitrarily bent monopole, above an infinite, perfectly 
conducting ground plane, in a Hallen-type integral which is founded 
entirely upon the familiar scalar potential and magnetic vector potential. 
These integral equations are then solved using the direct integration 



9 



technique. Throughout the thesis it is found that the coupled integral 
equations obtained by the Hallen approach can be readily obtained with- 
out employing an auxiliary scalar potential. How the various arbitrary 
constants, which arise in the development, are related to the electric 
scalar potential is clearly shown. 
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II. DERIVATION OF MATRIX EQUATION 



In this section the system of integral equations for the current 
distribution on an arbitrarily bent monopole, above an infinite, perfect- 
ly conducting ground plane, driven at its base by a delta gap generator 
is derived. From the equations, an operator is defined such that the 
method of moments can be applied to find an approximate representation 
of the operator by a matrix. The solution for the desired current dis- 
tribution is obtained by inverting the matrix. 

A. FORMULATION OF THE PROBLEM 

Figure 1 defines the coordinate system. The antenna consists of two 
thin wire segments with a single junction. The segments have radii a^ 
and a g and lengths h^ and h . The primary source is the field in the 
gap which is maintained at V^/A volts per meter by the generator. The 
induced source is the current on the dipole. The solution is carried 
out assuming the wires are perfect conductors. With this assumption the 
tangential component of the total electric field at the surface of each 
wire segment must vanish. Due to axial symmetry of each wire segment, the 
phi-component of each electric field must vanish. Hence for each wire 
segment the following is obtained: 




( 1 ) 



E S ( »S> + E 3 ( Ps> ' 0 



( 2 ) 



where E^(r g ) I s the tangential component of the induced electric field 
at the surface of the z-directed wire and E^(r g ) is the tangential 
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component of the primary electric field at the surface of the z-directed 
wire. Equation (2) is similarly obtained for the s-directed wire. 



B. DERIVATION OF THE COUPLED INTEGRAL EQUATIONS 



The Halldn equations are derived from a classical viewpoint for the 
case of a cylindrical antenna driven at its center by a so-called slice 
generator. [ Ref . 9, 18, 19]. The general procedure for developing a 
Hallen type equation is as follows: 

(1) Determine the magnetic vector potential component from the 
potential integral equation for A. 

(2) Determine the magnetic vector potential component along the 
cylinder from the inhomogeneous differential equation which the magnetic 
vector potential component satisfies. 

(3) Write the desired integral equation by equating the two ex- 
pressions obtained in steps one and two. 

(4) Apply appropriate boundary conditions on the surface of the ■ 
cylinder to the equation obtained in step three. 

1. Vector Potentials 

Consider the antenna shown in figure 1. Assuming harmonic time 
dependence (e^ Wt ), the vector potential due to the antenna current is 



A(r) 



V 

4tt 



J(r') 

J 



-jk | r-r ' | _ 

dv(r') 



(3) 



A(P ) = 



4n 



; 



J(p’) 



,-jk|p-p 

F-P ' I 



dv ( p ' ) 



( 4 ) 



Due to the geometry the vector potential has only a z-directed and s- 
directed component which can be written on the axis of each wire segment 
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as 



A z (z) ” 4» { I 2 (z ' ) 

z' 

V s > ‘rff v> 



e -3KR 



R 



dz 1 




(5) 

( 6 ) 



s 1 

where R = | r-r ' | and P = jp-p ! |. R and P were evaluated at r ’ and p* 

equal to the radius of the respective wire segments with (f> ! equal to zero. 

2 . Non-Ilomogeneous Differential Equations 

It is now seen that A (z) and A (s) can not be determined di- 

z v s 

rectly from equations (5) and (6) without knowledge of I(z') and I(s ! ). 

A differential equation is now derived which can be solved for A^(z) and 



A S (S) 



E = -jtHA + 



grad (div A) 



(7) 



ju)|j,e 

The electric fields produced by the currents are , according to equation 



(7), 

V (7) = jWA z (7)U z + grad [div(A z (7)uJ] 

E^CP) = -jo ) A s (p)U s + grad [div^ (p )uj] . 

Having applied the gradient operation in circular cylindrical coordi- 
nates and with axial symmetry eliminating the phi-component of E^ 1 (r) 

and E 1 (p) > the fact that the total tangential electric field on the 
s 

surface of .the antenna must vanish is applied. The following results 
are obtained: 



( 8 ) 

(9) 




A 

z 





E 




a r a 

3z L3y 




(10) 




A (p ) 
s s 



= iisi 

(JU 



E P (p 
s vr s 



> V p s>]- 



( 11 ) 
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where and denote any point on the z and s-directed wires respec- 
tively. Equations (10) and (11) are two coupled partial differential 
equations because of the last term on the right in each equation. A 
Hallen solution of these two coupled equations is needed. 

3. Solution of Coupled Differential Equations 
Equations (10) and (11) are of the form 



dx 

The general solution to equation (12) was contained in Ref. [23]. It 
was found to be 



00 



+ k f (x) - u(x) 4- 



dv (x) 
dx 



( 12 ) 



f (x) - A cos kx + B sin kx -f ~ 



px r 

, L u(x/) 






^ “q 



dV(x'V 
“dx' . 



sin k(x-x 7 )dx # . 

(13) 



The integral portion of equation (13) is zero at x equal q and its deriv- 
ative is zero at this position. One integration by parts results in 



px 



° k . -q 



^7 ~ ~ sin k(x-x / )dx / - -V (q) sin k(x-q) 



+ k 



V(x / ) cos k(x-x # )dx 7 



(14) 



°k =q 



which, when put into equation (13), gives 



f(x) = A cos kx 4- B sin kx - *t“ V (q) sin k(x-q) 




U(x') 

/= q 



sin k(x-x 7 )dx / -f 



px 



°x -q 



V(x 7 ) cos k(x-x , )dx / . 

(15) 



Since V(q) is fixed, the first three terms are combined and new constants 
obtained are 



a = A 4- — V (q) sin kq, and 
b = B — ~ V (q) cos kq. 



(16a) 

(16b) 
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Combining equation (16) with equation (15) gives 



ox 



f(x) - a cos kx + b sin kx -f 



U(x / ) sin k(x-x / )dx' 



~q 



ox 



V(x^) cos (x-x^dx' 



(17) 



~q 

With equation (17) available the solution to equations (10) and (11) is 
'determined with 

a ..2 

A y( r s ) = V( x )> V E z P ( r s ) = U ( x ) for eq uation ( 10 ) 
a -l 2 - 

and ~bw A w^ P s- ) = ^r E c P (P c ) - U(x) for equation (11) 



uu s xr s' 



which gives 



A (z) = a cos kz 4- b sin kz -f — 
z z z uu 



oz 



4n-z 



E^(n) sin k(z-n)dn 



* i; k v-> 

o 



cos k(z-n)dn and 



(18a) 



ik 

A (s) - a cos ks + b sin ks + - — E (n) sin k(s-n)dn 

S S S CD I S 

^ =s o 



" ^ L v n> ] cos k(s " n)dn * 



(18b) 



The v" A (z)J means A^ r (r) is differentiated with respect to y and then 



J3_ 

9yL y ' ' J y 

has to be evaluated along a z-directed line on which it is dependent. 

This is similarly true for — - j^A^(s)J. The arbitrary constants of equa- 
tion (18a) are a , b and z and any one of these three may be specified 
z z u 

directly, but the other two must comply with physical boundary condi- 
ditions. This applies to equation (18b) and constants a^ , b^ and s^ 
also. 



v 
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4. Scalar Potential 



With equations (18a) and (18b) one faces the problem of what to 
do about the constants of integration. One can enforce the boundary 
condition that the scalar potential must be continuous at each wire seg- 
ment junction to eliminate a constant or relate two constants. To obtain 
such an expression is the next task. 

Differentiating equation (18a) with respect to z and (18b) with 
respect to s according to Leibnitz 1 s rule [Ref. 24] gives 

a . , x , _ v 3 r jk 



3z z 



A (z) = -ka sin kz + kb cos kz + — E (n) sin k(z-n)dn} 

-7-7 -7 92 1 (JO J ^ * 



z 
n=z 



z 

f 3 



3y 



Ay(z) + k kjy Ay(n)] sin k(z-n) 



dn 



(19a) 



n=z 



and 



- A (s) = -ka sin ks + kb cos ks + — {— E (n) sin k(s-n)dn} 

<=■ ^ v o <= 9s L to J c J 



3s s 



J s 

n=s 



3w w 



s 

'Ife 



A (s) + k t— A (n)] sin k(s-n)dn. 



dW W 

n=s 



(19b) 



The Lorentz condition relates the magnetic vector potential A to the 
electric scalar potential and yields 

" • i 2 



- •Jk- A ( r} = — - A (r) 4- -r~ A (r) and 



3z z 



_ J* 



^ (r> Sy 

”- 2 a - 3 - 

M = aTT A T7^ P ^ + TZ • 

(0 oW W 3 S S 



(20a) 

(20b) 



Rearrangement of equations (19a) and (19b) gives 



- ^k ~ 4 > Z ( z ) = -ka sin kz + kb cos kz + -r-^- p— E (n) sin k(z-n)dn} 

2 >4 rr * f,\ -7 * 

z 

[[— 



dz 1 W j Z 

n=z 



+k |[— A^(n)] sin k(z-n)dn 



(21a) 



n=z 
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and 



2 

J k A S 



03 



4> b (s) = -ka sin ks + kb cos ks + — ^ E (n) sin k(s-n)dn} 

S S 3S 1 03 s J 



n=s 



+k 



A w (n)] sin lc(z-n)dn. 



J 

n=s 



(21b) 



which is valid along the z- and s-directed lines respectively, 
z 

<J> (z) is the total scalar potential along the z- directed wire 

s 

and <j> (s) is the total scalar potential along the s-directed wire, both 
of which are due to all charges on the wire segments, but the potentials 
are expressed in terms of the constants of integration. The partial 
derivatives expressed in equations (21a) and (21b) must now be evaluated. 

Going back to equation (13) it is seen that the integral portion 
of (13) is zero at x equal q and its derivative is zero at x equal q. 

Now in the derivation these forms are used 
x x 

~ J U(x^) sin k(x-x^)dx^ + j V(x^) cos k(x-x^)dx^ . (22) 



x =q 



x'=q 



Their relationship to the antenna problem is given below equation (17). 
Now applying these facts the result is 



_3 

3z 



{— E (n) sin k(z-n)dn} = 0 for all z 
1 (i) J z J o 



n-z 



{— E (n) sin k(s-n)dn} = 0 for all s 



(23a) 



(23b) 



n=s 



This result implies that one may choose s and z as desired, but at the 

o o 

same time doing this in a manner to lessen subsequent calculations. Once 
s^ and z^ are determined the other two constants of integration in each 



equation must meet boundary conditions. 
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C. BOUNDARY CONDITIONS 



The application of these necessary boundary conditions to the cur- 
rents on the thin wires should provide the additional equations to yield 
a unique solution for the current distributions. 

1 . At The Junction 

As has been stated by Taylor [Ref. A], there are essentially 
two conditions which must prevail at all junctions: 

(1) Kirchof^s current law must hold. 

(2) Continuity of the scalar potential must exist. 

The first of these requires that there be no buildup or accumulation of 
charge at the junction. Stated differently, the total charge entering 
the junction must be equal to the total charge leaving the junction. 

The second condition requires that the potential on each of the wires at 
the junction be the same. There can be no jump or discontinuity in sca- 
lar potential at a wire junction. 

2 . At The Free Ends 

According to previous considerations of thin wire theory [Ref. 1, 
2,4], the current must vanish at the free ends of all wires in the con- 
figuration. 
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III. SIMPLIFICATIONS ARISING FROM GEOMETRICAL SYMMETRY 



Any symmetries in the geometry of the antenna elements are reflected 

in the structure of the generalized impedance matrix Z . Judicious 

mn 

use of these symmetries can result in considerable reduction in both 
computing time and the number of memory locations required for a solution. 

A. IMAGE SYMMETRY 

Referring to figure 1 the antenna has image symmetry about the plane 
z = 0 through the midsection of the antenna. The boundary condition at 
a perfect electric conductor is zero tangential components of electric 
field. An element of source plus an "image** element of source, radiating 
in free space, produces zero tangential electric field over the plane 
through the midsection of the antenna. According to uniqueness concepts, 
the solution to this problem is also the solution for a current element 
adjacent to a plane conductor. The necessary orientation and excitation 
of image elements is summarized in many texts. The main point though 
is that symmetry must be maintained in the image problem. Once image 
theory is used and the resulting symmetry is recognized, all that is 
needed is to find the kernel for each "image" current element and com- 
bine it algebraically with its corresponding "non-image" current element 
kernel which reduces the number of required equations by a factor of two. 

B. TOEPLITZ PROPERTY 

The Toeplitz property of the Z matrix can be used to realize 

mn 

considerable economy of computing time. An element of the matrix depends 
only on the absolute value of the distance between the source point m 
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and the field point n if there is no bend involved. Since there are 
only N distinct distances involved, only N different quantities need 
be calculated for a matrix containing N x N elements. These N quanti- 
ties are calculated first and the computer is instructed to arrange 
them in the proper locations to form the complete matrix. 

The antenna has a bend so only certain portions of the matrix can 
be obtained jasing this property. It is the image of the bend that pre- 
vents the entire matrix from exhibiting this property. 
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IV. NUMERICAL SOLUTION 



One of the most general and most useful techniques for solving pro- 
blems in antenna theory is the direct integration technique. A detailed 
discussion of the technique has been presented by Harrington [Ref. 16]. 
Since this basic technique is employed throughout the remainder of this 
thesis, a short quantitative discussion using the notation of Harrington 
is included at this time. 

A. METHOD OF MOMENTS 

The equation to be solved is of the form 



where L is a linear operator, g is the tangential component of the applied 
field on each conductor and f is the unknown current to be determined. 
Using the terminology of Harrington [Ref. 17] it is evident from the 
uniqueness theorem that the problem is deterministic. Since L and g are 
given, it is a problem of analysis. 

A unifying concept for the solution of such problems is afforded 
by the method of moments. In order to apply the method, it is necessary 
to define: 



domain of L. These known functions must be linearly independent and 
chosen such that a reasonable approximation to f can be obtained by a 
finite series expansion of the form 



L(f ) = g 



(24) 



(1) A set of expansion or basis functions f^, f^ 



f in the 
n 



N 




( 25 ) 
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where the a ! s are unknown constants to be determined. If equation (25) 
n 

is substituted into equation (24), one obtains 
N 

l a n L(f n } = g ‘ (26) 

n=l 

(2) A set of weighting or testing functions W^, ^in+l 5 

. . . ,W^ in the range of L and linearly independent must be chosen. If 
the same set of functions is used for expansion and testing, the pro- 
cedure is known as Galerkin’s method and is equivalent to the Rayleigh- 
Ritz variational procedure [Ref. 17, 25]. 

(3) An inner product suitable for the problem. For most electro- 
magnetic field problems, a suitable inner product is 

f " 2 

<W,f> = f (n)W(n)dn* (27) 

n l 



where n^ and n^ are chosen over the region of interest. The resulting 
expression is 

N 

I a <W ,L(f )> = <W.g> , m = . (28) 

I Li ill 11 Lil 



Inspection of equation (28) shows that it can be written as a matrix 
equation by defining the following matrices 



U ] = 

mn 



<W 1 ,Lf 1 > <W 1 ,Lf 2 > ... <W 1 ,Lf N > 

<W 2 ,Lf x > <W 2 ,Lf 2 > ... <W 2 ,Lf N > 



<W H ,Lf l > <W N ,Lf 2 > ••• < V Lf N > 



(29) 



= 


a l' 


, and [V ] = 
m 


\ 

A 

OQ 

V 




a 2 




<W 2 ,g> 




• 

• 

. V 







( 30 ) 
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The matrix equation may now be written as 



U J [a] = [V 1 

mn n m 



( 31 ) 



By analogy with multiport network theory, [z ] is called a generalized 

mn 

impedance matrix and [V ] a generalized voltage matrix. Only if [z ] 

in mn 

is nonsingular does its inverse exist and the desired solution is 



[a ] = U ] 1 [V ] . 

n mn m 



(32) 



Since the a f s are then known, the solution for the unknown current f 
n 

is determined by equation (25). 

The expediency of obtaining a solution for f depends not only on L 

and g, which are usually determined by the physical problem, but also 

on the functions f and W . Extreme care should be taken in choosing 

n m 

these two sets of functions. Experience with problems of the general 
type considered in this thesis suggest that only certain particular 
choices for the expansion functions and for the testing functions yield 
equations which are amenable to solution. In particular, when the test- 
ing functions are chosen to be Dirac delta functions, the resulting 
technique is usually referred to as point matching. 



B. EXPANSION FUNCTIONS FOR THE UNKNOWN CURRENT 

There are several popular choices for the expansion functions which 
have been considered in the literature. Some of the more common of these 
are piecewise constant, piecewise linear, piecewise parabolic and piece- 
wise sinusoidal [Refs. 20, 21]. All of the aforementioned fall under the 
general label of the method of subsections. This method has the antenna 
contour divided into a number of intervals or subsections where the 
boundary condition of equation (1) and (2) is to be applied. The lengths 
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of the intervals need not all be equal. For certain types of problems 
it could be advantageous to weight the lengths of the intervals. In order 
to keep the formulation as simple as possible, however, the present 
solution uses N-l intervals of equal lengths and two intervals of half 
the length of the N-l intervals. These half subsections were to allow 
a match point to be placed at the ends of the antenna and one on the 
bend of the antenna. In addition to these piecewise functions is another 
family of expansion functions which are commonly called entire-domain 
expansions. A discussion of the use of these entire-domain functions 
will not be attempted in this thesis. 

All of the calculations and results presented in this thesis are 
obtained by the method of subsections in conjunction with point match- 
ing. In particular a series of piecewise constant and piecewise sinu- 
soids are used to represent the unknown currents. The resulting 
equations are forced to be satisfied at a set of discrete points called 
match points. Using Dirac delta functions as weighting functions ef- 
fectively requires that equation (24) be satisfied at N points on the 
contour of the antenna. 

Harrington and Mautz [Ref. 2] have solved the single straight wire 
problem using three different procedures, all within the scope of the 
method of moments. In order of increasing computational complexity, 
these include: 

(1) Point matching with pulses as basis functions 

(2) Point matching with piecewise linear basis functions 

(3) Galerkin's procedure with piecewise linear basis functions. 

As the number of subsections increases, the second solution converges 
about twice as fast as the first. With antenna segments less than A/10 



25 



in length, no significant difference in results is observed between 
Galerkin f s procedure and the second procedure. Using piecewise sinu- 
soids is equivalent to the second method in speed of convergance and 
and accuracy, but piecewise sinusoids do not require the segment length 
to be less than A/10 in length to obtain this accuracy and speed. 

Using Dirac delta functions as weighting functions gives 



W = 6(z - z )U 
m m z 

A 

W = 6(s - s )U 
s ms 



(33a) 

(33b) 



in which is a unit vector directed along the positive z-axis and U g 
is a unit vector directed along the positive s-axis. Using the piece- 
wise constant subsection gives the following expressions for current 



N 

i( Z ) = y i p ( z )u 

L . n n z 
n=l 



(34a) 



I(s) -VIP (s)U 
q q s 
q=l 1 n 



(34b) 



where 1^ and 1^ are the complex constants to be determined and the pulse 
functions are constant over the increments. This statement in equation 
form is 

1 z <z<z 



P (z) = 
n 



n n+1 
0 elsewhere 



(35a) 



p q (s> - 



1 s <s<s . i 

q q+i 



(35b) 



0 elsewhere 

Similarly, the piecewise sinusoid expression gave the following current 
equations 

N I sin k(z f ,.-z f ) + I sin k(z f -z T ) 

I(z) = y 2+i 2+1 2 — p ( Z )U (36a) 

L n z 



n=l 



sin k(z ' . -t — z 1 ) 
n+1 n 
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Q I sin k(s’ -s 1 ) + I - sin k(s ! -s ! ) 

Ks) - l -3 3±i a±l 2_ P <s)u (36b) 



q = i 



sin k(s 1 . -s 1 ) 
q+1 q 



C. LOCATION OF MATCH POINTS 

There has been little quantitative study done on the precise effect 
of the location of the match points on the solution; however, it has 
been the experience of this author that the location of these match 
points significantly affects the accuracy and convergence of the solution. 
The number of match points and their location depends on such factors as 
the form of the integral equations (Pockington or Hallen) , the type of 
current expansion used (pulse, sinusoids) and the allowable time and 
core for computation. 

Extensive calculations have been carried out [Ref. 22] where the 
match points are located at ends of each current subsection. The results 
of these calculations indicate the solutions are quite stable and reason- 
ably accurate even for problems where a large number of match points is 
used. 

Other investigations have indicated that a good choice for the loca- 
tion of the match points is at the center of a current subsection. As 
the author discovered these conclusions have to be tempered with other 
information. One question is how and what kind of experimental verifi- 
cation is used. If a scattering cross section is used, the demonstrated 
results may differ from those obtained for an excited antenna. The reason 
is the cross section calculations effectively average the source and 
current distributions on the structure. Thus, a convergence study 
pertaining to input admittance, using scattered results, may be somewhat 
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inconclusive since the segmentation and match point dependence are 
directly related to source model and numerical convergence. 

All of the results presented in this thesis are the result of locating 
the match points at the center of a current subsection for piecewise con- 
stant expansion and at the ends of the current subsections for piecewise 
sinusoid expansion. 
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V. APPLICATION OF THIN WIRE THEORY TO AN 



ARBITRARILY BENT WIRE MONOPOLE 



A. GENERAL CONSIDERATIONS 

A general theory for determining the currents induced on an arbitrar- 
ily bent wire monopole has been outlined. The purpose of this section is 
to illustrate how to obtain the exact solution. The configuration to be 
studied is composed of two straight nonorthogonal elements intersecting 
at a single junction which are driven above a perfect infinite ground 
plane. The problem is developed in a manner which allows general var- 
iations in many of the fundamental parameters. For example the lengths 
of the individual elements may vary, the angle between the individual 
elements may vary and the radii of the elements may vary. It should be 
noted that all variations in parameters are still subject to the basic 
assumptions and restrictions imposed by thin wire theory. These restric- 
tions do not greatly inhibit the usefulness of the investigation since 
much general knowledge can be obtained from the thin wire approach. 

B. COUPLED INTEGRAL EQUATIONS 

The coupled integral equations for magnetic potential which were 
developed in a previous section are 



z 



A (z) = a cos kz + b sin kz + - 
z z z 




E (n) sin k(z-n)dn 
z 



n=z 



z 



o 




(37a) 



n=z 



o 
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i k 

A (s) -a cos ks + b sin ks + - — 
s s s 0 ) 



E (n) sin k(s~n)dn 
s 



n=s 



" r { ^(n)]} coS k(s “ n)dn 



(37b) 



n=s 



and the scalar potential equations are 



M 2 



- <|> Z (z) = ~ka^ sin kz + kb^ cos kz + k [ {— [A^(n)]} sin k(z-n)dn (38a) 



n=z 



and 



ik 2 z 

— d> (z) - -ka sin ks 4- kb cos ks + k 

03 T V S S 



{-— [A^(n)]} sin k(s-n)dn. (38b) 



n=s 



The constants a^, b^, a^ , and b^ must assume values which render equa- 
tions (37) and (38) in compliance with the following boundary conditions: 

(1) The scalar potential, <j> , must be continuous along the wire 
segments except at the delta gap generator and in particular the scalar 
potential must be continuous at the wire junction. 

(2) The magnetic potentials, A g and A^ , must be continuous along 
the wire segments and at the wire junction, 

(3) Current must vanish at the free end of the wire. 

(4) The axial currents must satisfy Kirchoff f s current law at 
the wire junction. 

The magnetic potentials, A^ and A^, developed earlier, are. 



and 



A (z) = “t— 
z 4 tt 



A (s) = 
s 4 tt 



z ' 

r 



“jkR 

V^V-* 2 ’ 



- jkP 

I s ( s ')^-p— ds* 



(39a) 



(39b) 
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where R = | r - r 1 | and P = |p - p 1 | . R and P were evaluated at r 1 and 
p’ equal to the radius of the respective wire segments and cj) 1 equal to 
zero. Using these values for |r - r* | and |p — P 1 | and assuming the sur- 
face current flows as an equivalent line current along the s and z axes 
is equivalent to using the reduced kernel. This will be proven in Ap- 
pendix A along with conditions when the approximation is no longer valid. 

Looking back at the development of equation (37) it is seen that the 
mathematical formulation was given without any physical insight as to 
where in space the equation was valid. In going from equations (10) 
and (11) to equation (18) there is a built-in assumption not expressly 
stated. This assumption is that the radius vector is fixed at some value 
which results in equation (37) being valid only along a line parallel to 
z and s directed wires. This result will be fully shown in Appendix B. 
Since equations (37) are valid along any line parallel to the z and s 
directed segments, one can require them to be along the surface of the 
wire segments which is along a line (a,o,z) and (a,o,s) respectively. 

This surface is also where the tangential E field is zero. Since the 
source points and field points could not lie along the same line, the 
reduced kernel has to be used with its equivalent line current flowing 
along the axis of the two wire segments. 

Figure 1 defined the co-ordinate system which is composed of two 
co-ordinate subsystems. They are the (X,Y,Z) and (P,W,S) co-ordinate 
systems. The kernels for the wire segments are then expressed in these 
co-ordinate systems . 



-jkR 

K (x,y,z,z') = - — - — 
Z R Z 

-jkR S 

K s ( x ,y»z» s ') = - — - — 



K*(x,y,z,z’) 

Kg(x,y,z,sp 





(40a) 

(40b) 
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(40c) 



G z (p,w,s,z') = 



G (p ,w,s ,s 1 ) = 



-jkr 



-jkr- 



-jkr Z 

x , , x e x 

G z (p,w,s,z.) - “ — 

r . 
x 

« 1 S 

-ikr, 

~X. IN e J X 

G s (p,w,s,s.) = — 

r . 
x 



(40d) 



R = {(x - x') 2 + (y - y') z + (z - z') z } 



• \ ^ 



. t \ 2i 2 



r = { (p - p')^ + (w - w') 2 + (s - s'/'}^ , in general. 

In order to account for the images due to a perfect ground plane one 
must include the image current kernels which are denoted above by the 
superscript i. 

With reference to figure 2, 



i n ^ \ 2 



and 



0 



A(r) = U { I (z f )K (x,y,z,z T )dz ! + 1 (z ! )K 1 (x ,y , z, z 1 ) dz ! } 

z 1 / z z zxz xx J 



h 

s 



+ U 



z ! =-h 

.O 1 2 



s’=0 



I (s')K (x,y ,z ,s ' ) ds 1 - U . 

s s’ si J 



I (s')K^(x,y,z,s')ds' (41) 



sis 

=-h 
1 s 

0 



i i 



A(p) = U { I (z')G (p,w,s ,z')dz' + 



h 

r s 



z =0 



I (z! )G 1 (p,w,s,z!)dz.' } 
z l z l i > 

z ! =-h- 



r 0 



+ U s J I s ( s ') G s (P> w > s > s ') ds ' “ u 

s ' =0 



si 



I (s '. ) G 1 (p ,w s ,s ! )ds ! . (42) 

S X S X X 



s!=-h 
x s 



When the U . and U unit vectors are broken into their respective U 
sx s z 

and U components, they give the required A (r) and A (r) . Similarly 
y z y 

/V /V A A 

breaking and into their respective and components gives 

A g (p) and A^(p) . These four magnetic potentials are needed in equation 



37. These four required magnetic potentials are: 
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A (r) 

V z 



h 

s 



I (z')k (x,y , z , z ' ) dz ' + cos0 



I s (s')k‘ (x,y,z,s ' )ds ' , (43a) 



z'=0 



s ' =0 



— — A (r) 
y y 



A <p) 
y s 



= sin 0 j 1^ (s ' )k^(x,y , z, s ' ) ds ' , 

s '=0 
h 

( s 

= I 1^ (s ' ){ G^ (p ,w, s , s ' ) ds ' + cos 0 G^(p,w, s ,s ' )} ds ' + 
s ' =0 

h 



(43b) 



cos 0 | 1^ (z ' ) G Z (p ,w, s , z ' )dz' 



and 
4 TT 



A (p) 

p w 



z'=0 

= -sin 20 



h 

s 



I (s ' )G (p ,w,s ,s ' )ds 1 - 



sin 0 



s '=0 
h 

f Z z 

I (z*)G (p,w,s , z' )dz' . 
z'=0 



(43c) 



(43d) 



Placing equations (43) into equations (37) gives 



11 11 
r z r S 

~ (z T )K Z (x ,y , z, z ! ) dz T + cos 0 (s T )K S (x,y , z,s T )ds f = a^cos kz + 



z T =0 



s ! =0 
z 



b^sin kz + j E z (n)sin k(z-n)dn - — j {y~[A^(n)]} cos k(z-n)dn (44a) 



and 

h 



n~z 



n=z 



^ 1 

I (s T ){G (p,w,s,s f ) + cos G (p ,w,s ,s T ) }ds T + cos0 



h 

r z 



I z (z T )G (p ,w,s , z T ) dz 1 



s T =0 



z T =0 



jk 



a cos ks + b sin ks + 
s s co 



i s 

E g (n)sin k(s~n)dn - — j f j^tA^n)] } cos k(s-n)dn. 



n=s 



n=s 



(44b) 



Examination of equations (44) indicated that the partial derivatives 

of A (r) and A ( p) are needed. The results for A (r) will be shown and 
y w K y 



A (p) will follow in a similar manner, 
w 
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The form of A^(r) is 



A y (r) = 



,-jkR 

~R 



and 



9A (r) 9A (r) 9R 

y = __y 

9y 9R 9y 

9A (r) 



9y 



_ (1 + j k R) (y - 



,* -JkR 

y )e ■’ 



R“ 



(45a) 



(45b) 



Similarly, for A (p ) , 
w 



9A (p ) 
w 

9w 



(1 + jkr)(w - w')e 
3 



(46) 



r 

Placing equations (45) and (46) into equations (44) and examining 
the results closely, it is seen that the double integrals have an inner 
integral of the form 



where 



[ + -y] e J cos k(3-n)dn 

J R R 

1 

R = {a 2 + g(a) + [n - f(a)] 2 } 2 



(47) 



and f and g are any function of the parameter a. Reference [23] gave 
the exact solution for this integral which reduced the double integral 
to a single integral and thereby greatly reducing the computer inte- 
gration time and numerical errors. The exact derivation will be shown 
in Appendix C. Now place the result of Ref. [23] into equation (44a). 
With rearrangement of equation (44a) the following equation is obtained: 
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[— ] a cos kz + [— ] b sin kz + [ - 71 J — ] I E (n) sin k(z - n)dn - 
\j z y z yu | z 



n=z 



4tt 



x / IVY/ , . r [n-f (s' ) ]cos k(z-n) - jRsin lc(z-n)i % 

I (s') (K- X (x,y,n > s'){ 1 o 2 ) J ds 

^ J , ! R(R - [n-f(s')n 

s =s n=z 

h h 

Z r S 

I (z 1 )K Z (x>y , z 5 z f ) dz f - cos 0 [ I (s 1 ) K S (x ,y , z y s 1 ) ds 1 = 0 
z I s 



( 48 ) 



Similarly , for equation (44b): 



[— ] a cos ks + [— ] b sin ks + [AlLlli] 
y s y s ya) 



E (n) sin k(s - n)dn + 
s 



n=s 



4rsin 20 
V 



4rsin 0 



) ds ' + 

n=s 



f I (s') (c i (p,w,n,s')( |n - f(s ' )ll: ° s , k(s - n) - lrSkn k(s - n ) | 

l ,J r(r - [n-f (s') p) 

s s 

h s° s - 

I (s') (o z (p , w , n , z ’ ) l-is.- - »* ■' -) j£°S Sk fe ” ) . } ] ds . _ 

, S r(r - [ -f(z')] ) 

s '=s n=s 



I (s') [G (p,w,s,s ') - cos G (p,w ) s > s')]ds' - 
o s s 



cos 0 I I z (z f )G (pjWjSjZ^dz 1 = 0 



(49) 



For a structure driven by a delta gap source of emf, V, on an element 
at a location displaced g units above a ground plane the following equa- 
tions are obtained [Ref. 20]: 

z 



V (z) = & 

Z 0) 



/ 

-VS(n-g) sin k(z-n)dn 
n=z 



v z (z) 



= _ jjs. 



(0 



[V sin k(z-g) U(z-g) ] . 



(50a) 



(50b) 
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It is desirable to obtain a representation of the emf when the 
structure is base driven, i.e., when g is zero. A structure above a 
ground plane, regardless of the nature of its excitation, always gave 
rise to zero scalar potential everywhere on the ground plane. This is 
a direct result of the uniqueness theorem and boundary conditions at a 
perfect conducting surface. Thus, it is necessary that, if the element 
and its source are to be equivalent, the impressed emf must be sym- 
metric; otherwise, the currents and charges in a free-space ’’equivalent" 
of the structure above a ground would not satisfy imagy theory require- 
ments . 

It is from this discussion that the modification to equation (50b) 
is found. When a base-driven structure is being analyzed, one sets g 
equal to zero and V equal to V/2. If V/2 were not used in place of V, 
one would be treating an antenna driven by 2V when the free-space 
"equivalent" structure is formed. These results are 

V (z) = - — — sin kz . (51) 

z to 2 

Since the s-directed segment is not driven by a separate generator, this 
gives 

s 

V (s) = 0 or — [ [0] sin k(s - n)dn = 0 . (52) 

s 0) J 

n=s 

o 

Now the boundary conditions are applied. Due to symmetry of the 
charges on the structure and its image, the scalar potential must be 
zero at the point where the vertical element enters the ground plane. 
This implies 

. ,2 

- -L- <j) Z (0) = 0 (53a) 

which gives = 



0 



(53b) 



Next by condition one the scalar potential must be continuous at the 
wire junction which gives 

<j> Z (h z ) = <j> S (0) (54a) 

and results in 



-ka sin kh - kb + k 
z z s 



{— [Ay(n) ] }sin k(z-n)dn - — ~ cos kh^ = 0 (54b) 



n=z 



Conditions three and four require 



I Q = 0 and 



(55) 



i z = i s 
N 1 ' 



(56) 



Clearing all constants and applying all boundary conditions, the fol- 
lowing equations for I (z) and I (s) are obtained: 

z s 



a z 

[— ]a cos kz - f I (s 1 ) [K y U.y ,n ,s 1 ) (.Lq_ ^ ' l l kosJiU - n)_ - j E sin k (z- n) j ) Js , 

“ 2 i'.o R < R - „-0 



s ' =0 
h z 

- ^ sin kz - j 1^ (z* )K Z (x,y , z , z' ) dz ' - cos 0 
z'=0 



h 

s 



I s (s')K (x,y , z , s ' ) ds 1 



s ’=0 



= 0 



(57a) 



f-^-1 a cos ks + 
V z 



; ' ) f [n-f(s , )]cos k(s-n) - jrsin k(s-n) i I -j dg , 



r(r 2 - [n-f(s')] 2 ) 



I s (s’)(G 1 (p,w,n,s'){ 

s '=0 A n=0 

h s 

z, , x j [n-f (z 1 ) ] cos k(s-n) - jrsin k(s-n) 



+ [~M’]b s Sin ks + | I z ( z ') (G Z (p,w,n,z'){- 

z ' =0 



r(r 2 - [n-f(z')] 2 ) 



) ) dz ' 

n=0 



h h 

s ( s i 

I (s*)G (p,w,s,s f )ds f - cos 0 I (s T )G (p ,w, s , s ! ) ds 1 
S S J s s 

s ' =0 s ’ =0 

h 

z 

I z (z' )G Z (p,w,s,z' )dz' = 0 

z ' =0 



- cos 0 



(57b) 
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and 



[— ]a sin kh + 



y z z 







4tt 




0 



(57c) 



y s 



z 




0 



(57d) 





(57e) 



There are five unknowns a , b , a , I (z ! ), and I (s T ) to be found 

s s z z s 



and five equations which relate these quantities; hence, there is a unique 
solution to this system of equations for the five unknowns. [Ref. 20], 
Obtaining the current distributions, all the other desired quantities 
such as electric field, input impedance can be determined. 

C. NUMERICAL SOLUTION 

The system of integral equations given in the previous section must 
be solved in order to determine the unknown current distributions. The 
direct integration technique as outlined earlier is used to effect a nu- 
merical (discrete) solution for the current distributions. 

The first step in the direct integration technique is the choice of 
the expansion functions and the testing functions. As pointed out in an 
earlier section, the method of subsections in conjunction with point 
matching is a very convenient approach. The testing functions are Dirac 
delta functions since point matching is used. 

The expansion functions for the unknown currents are chosen to be 
piecewise constant for its ease of use and piecewise sinusoids for their 
improvement in convergence and required computation time. Symbolically, 
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the piecewise constant function is represented by 



N 



V 2 ’) - £, I n 1U(z ' 

n=l 


- z') - U (z ' 
n 


- z ; + i> ! 


(58a) 


i (s’) = y i s [u(s ' 

S q=l q 


- s') - U(s ' 

q 


- S q + 1> 1 


(58b) 



and where U(a) is the familiar unit step function 



U(ot) = 



! 1 a > 0 



0 a < 0 



(59) 



Consistent with the above description, the current can be written as 

N-l 

I (z') = I Z U(z' - z') + l I z [u(z'-z') - U(z'-z' )] + I z u(z'-z') (60a) 

z 1 Z o n n n+1 N N 

n=2 

and 

Q-l 

I (s’) = I z u(s! - s’) + l I®[U(s’-s') - U(s ' -s' ) ] + I®U(s’-s') . (60b) 

s 12 ^2 q q q+1 Q Q 



Similarly, one can write for the piecewise sinusoids the following ex- 
pressions : 

N I Z sin k(z' ,-z') + I Z ..sin k(z'-z') 

I (z') = 1-* ^ — {u(z ' -z ' ) - U(z'-z’ .)} 

z _ . « f • i \ n ni* 



n=l 



sin k(z ' , ,-z' ) 
n+1 n 



Q I S sin k(s' -s') + I s sin k(s'-s') 

Us’> - l -3 S±1 3±1 Mu(s’-s’) - n(s’-s' +1 )} 



(61a) 



q-i 



sin k(s ' , ,-s ' ) 
q+1 q 



(61b) 

With the assumption that the segments were all equally spaced, the term 

z'-z' is a constant. This allows a term delta (A) to be defined as 
n+1 n 



A = — 7 



sin k (constant) 



(62) 



AO 



II 



With this definition and equations (61) the current can be written as 



I (z) = 
z 



sin Y^(z^-z' )U(z^-z 1 ) + 



N+l 



sin k(z f -z T )U(z 1 -z') + 
n N 



(63) 






I Z {sin k(z ! -z f ) fu(z ! -z f )-U(z*-z')1+ sin k(z T -z’ ) (u(z f -z f )-U(z ! -z f )1 } 
> n n-1 v n-1 n n+1 v n n+l J J 



and similarly for the s current. With these expansion functions, the cur- 
rent distribution appeared graphically as shown in Figures 3-9. 
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Figure 3. Real and imaginary parts of the current distributions for the 
radius = .001588A, theta =0°, 36 segments. 
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Figure 4. Real and imaginary parts of the current distributions for the 
radius = .004763A, theta = 0°, 20 segments. 
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Figure 5 
radius = 





. Real and imaginary parts 
.001588X, theta = 0°, 30°, 



of the current 
45°, 60°, 85°, 



distributions for 
90°, 36 segments. 



the 
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Figure 6 
radius = 





Real and imaginary 
. 003175X , theta = 0° 



parts 

30°, 



of the current 
45° 5 60°, 85°, 



distribution for 
90° , 32 segments 



the 



45 



Figure 7 
radius = 





• Real and imaginary parts of the current distribution for the 
.004763A, theta = 0°, 20°, 45°, 60°, 85°, 90°, 25 segments. 
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Figure 8 
radius = 





• Real and imaginary parts of the current distributions for the 
.001588A , .003175A, .004763A, theta = 45°. 



4 7 



Length ( > Lergth ( X ) 



Figure 9 
radius = 





. Real and imaginary parts of the current distributions for the 
.001588A, .003175A, .004763A, theta = 90°. 
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VI. RESULTS 



The system of integral equations and the proposed solution techniques 
are quite general, and they allow variations in many of the parameters in 
the basic configuration. However, certain parameter variations appear to 
be more interesting than others. This study concentrates its efforts in 
three main areas: 

1. The ability of the numerical solution to give the result obtained 

by King in Ref. 26 when the angle between the two wires is zero. 

2. Current as a function of the angle between the two wires. 

3. Current as a function of w ire diameter. 

The antenna is assumed to be driven by a delta gap generator. The 

current distributions are presented in terms of a dipole vice a monopole. 
The real and imaginary parts of the current are plotted for various bends 
of the antenna. The total length is A/4 for all the figures with the an- 
tenna bent at A/8 to give two equal lengths. Although these results are 
the only ones presented, the computer program allows the antenna to be 
bent at any point and the antenna to have any length. 

Figures 3 and 4 show the current for two different radii at an angle 
of zero degrees for a comparison against King’s frill feed in Ref. 26. 
These figures also show the result from Ref. 27 and 28 against King. 

References 27 by D.C. Kuo and B.J. Strait at Syracuse University uses 
a modified Pockington formulation with thin-wire approximations. 4> is 
used instead of A and derivatives are approximated by finite differences 
over the same intervals used for integration. Galerkin’s procedure is 
used with triangle expansion functions which are non-zero only over four 
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consecutive segments. Successive triangles overlap every two segments 
except on the ends of open wires. Junctions are treated as open wires 
each having a two segment overlap at the junction. An extended operator 
is used because the first order derivatives of the triangle functions 
do not exist at certain points. 

Reference 28 by MB Associates and called the Antenna Modelling Pro- 
gram uses a Pockington formulation with thin-wire approximations. Point 
matching is used with the current being expanded as a sine, cosine and 
constant with an interpolation procedure so only N linear equations need 
to be solved. The extrapolated current for a given segment is forced to 
match the center current values in two adjacent segments. This interpol- 
ation procedure is also applied to wire junctions. 

Figures 5,6 and 7 give the currents for antennas of three different 
radii having 0, 30, 45, 60 and 90 degree bends. King’s input values for 
a 0 degree bend are also plotted. Figures 8 and 9 show the currents for 
different radii antennas plotted on the same graph for bends of 45 and 90 
degrees , respectively. 

Several general comments concerning the data presented may prove help- 
ful to the reader. All the results presented here are for the case of a 
one volt drive and assumes an open-ended perfectly conducting hollow 
cylindrical tube as the antenna. Extensive calculations have been carried 
out for bends of 0 to 90 degrees for each antenna radius given, but the 
plotting of all these angles on a single graph would obscure the overall 
view. It is believed that the data given is sufficient to give the reader 
an indication of the general variation of the current with respect to the 
parameter of interest. 
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VII, CONCLUSIONS 



In this thesis a new method is presented for analysis of an arbitrarily 
bent monopole. The arbitrarily bent monopole is a good choice because it 
contains cross- coupling terms in the Hallen formulation which are not 
present in the 0 or 90 degree case. The currents had to be obtained 
using the appropriate boundary value problem based on Maxwell’s equations 
which are formulated as a operator equation. The general procedure, 

"method of moments", is used to represent the operator by a matrix. Geo- 
metrical symmetries are utilized to simplify the matrix inversion. 

Numerical results are compared with data obtained by different theo- 
retical techniques and with experimental results where possible. Agree- 
ment is generally very good and this is interpreted as an indication of 
the accuracy of the solution. 

Certain parameters should now be concentrated on if the conclusions 
are to be used properly. These efforts should be in these areas: 

1. Choice of bases 

2. Antenna drive 

3. Discontinuity at the bend. 

As can be seen by equation (36) the sinusoid current expansion 
equations are continuous but their derivatives are not. This implies 
mathematically that there are pockets of charge at the end of each 
subinterval which physically are not there. If one is interested in 
the near field, these mathematical charges cause discontinuities in the 
near field. One solution has been to use a least square fit to the cur- 
rent values obtained. 
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The delta gap drive was chosen because it was an ideal model. Even 
though the delta gap gives a physically meaningless infinite input sus- 
ceptance associated with the knife edges, the answer is close to that 
achieved by experiment. The reason is that any method that approximates 
the current by a series of continuous functions will not reproduce the 
extremely sharp rise in amplitude near z = z 1 unless a very large number 
of terms is used. A better drive would be the frill feed as given in 
reference [29] which corresponds closely to the actual drive used on 
physical antennas . 

Although not plotted on the figures given, a discontinuity in cur- 
rent was obtained within one subinterval either side of the bend. This 
was to be somewhat expected because the formulation considers two straight 
wires joined at a point which gives a sharp bend where the cross section 
is abruptly changed. When one uses the approximate kernel, I(z) must be 
a slowly varying function in the vicinity of z = z 1 for the approximate 
kernel to be within an order of the square of the radius of the exact 
kernel. This is not true at the bend. It is this slow variation 
or fast variation that causes the various methods of iteration (Hallen, 
Gray , King-Middle ton) to require more or fewer terms in representing the 
current. 
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APPENDIX A: REDUCED KERNEL 



Equation (39) in the thesis must be investigated and R and P deter- 
mined in a more precise manner. A detailed development of the vector 
potential must now be presented. 

The vector potential due to a volume current density can be written 
as 



Looking at Figure 10, it is seen that the vector r locates the point at 
which A(r) is to be determined and r* locates the source point. With a 
symmetrically driven antenna whose wall thickness approaches zero, the 
vector potential has only a z component which can be written as 




( 1 ) 



2tt a 



• i i « i 




or integrating over r 1 gives 




d^'dz 1 . 



(3) 



Now due to symmetry 



I(z ! ) = 2TraJ z (a,z ! ) 



(4a) 



or 




(4b) 



Substituting (4b) into equation (3) gives 




( 5 ) 
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Figure 10. Antenna Geometry assuming a surface current in 
the z direction. 
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where r 1 is constrained to remain on the surface of the cylinder. 

Appealing to basic vector algebra and trigonometry and being aided 
by Figure 10 one can write r and r ! in the following form: 

r = xU 4- yU + zU = r cos <j>U + r sin <j>U 4- zU (6a) 

x yz x y z 

AAA A A A 

r 1 - x f U 4* y 1 U +z ! U = r'cos d> f U 4- r 1 sin <f> f U 4- z ! U . (6b) 

x y z x y z 

Now |r - r 1 | is written as 



1 

|r - r'| = [(x - x’) 2 + (y - y’) 2 + (z - z’) 2 ] 2 

1 

= [r 2 + (r') 2 - 2rr'cos(<}> - <j>’) + (z - z') 2 ] 2 (7) 



If |r - r ! | is evaluated at r 1 equal to the radius of the z-directed 
wire and <j> equal to zero then 



|r - r f | = [r^ 4- a^ - 2ra cos <f> ! + (z - z 1 )^]^ 

r “a 

<f) -0 

Trigonometry gives 

2 

1 - cos 23 = 2 sin 3 



( 8 ) 



(9) 



Applying equation (9) to (8) with r = a gives 

1 

| r - r ' | = [4a 2 sin 2 ^) + (z - z') 2 ] 2 . (10) 



Going back to equation (5) and inserting equation (10) gives 

2tt r/ 2 . 2,<f>') ^ , ,.2.1/2 

r -jk[4a sm (-=— + (z - z ! ) ] 



A z <a - Z) ■ 47 






2ir 



. , = 0 [4a 2 sin 2 (|~) + (z - z ' ) “ ] 






One sees immediately that at z equal z ! and <p' equal- 0 or 2 tt the inner 
integrand is singular; however, the integrand is logarithmically singular 
and is integrable. Even though it is integrable, evaluation of the in- 
tegrand by numerical techniques is extremely difficult. This is one 
major reason why a reduced, average or approximate |r - r ! | is found. 
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For | z — z T | >> a and for thin wires where a << A the following 



approximation can be obtained. 



r - r'| = [(a - a ') 2 + ( 2 a sin l^]) 2 } 1 ^ 2 = |z 



( 12 ) 



Even this for thin wires is not a good approximation where z = z 1 . Ob- 
serving that 2 a sin (■——) is the distance from a point (a, 0 ,z f ) to a 
point (a,<J>*,z), it is seen that both points are on the same circle of 
radius a in a constant z ! plane. Looking at Figure 10 it is seen that 
2 a sin (y-) was a distance that assumed values from 0 to 2 a as (J ) 1 
varied from 0 to 2i t with a mean value of a. Placing the mean value into 
equation ( 11 ) and integrating gives 



A (a, z) = — 
z 4 tt 



I(z' )- 



r jk[a 2 + (z - z ’) 2 ] 1/2 
[a 2 + (z - z ') 2 ] 1/2 



dz 1 



(13) 



It is restated that equation (13) is valid only for a << A and 
a << h. Now going back and using Figure 11, equation (5), and assuming 
the current to be a filament along the z axis, the vector potential 



A^(r) at a point (a,0,z) would be 



A z (a,°,z) - 



I( z'Y 



r 2 . , ,.2 1/2 

~3Ma + (z - z') ] 

r 2 4 . / , x 2 , 1/2 

[a + (z - z') ] 



dz' 



(14) 



which is exactly the same equation as equation (13). This is a most 
important point. Using the reduced kernel is the same as treating sur- 
face current on the cylinder as if it were an equivalent line current 
along the cylinder axis. 
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Figure 1 1 . 
along the 



l 




Antenna Geometry assuming a filamentary current 
axis. 
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APPENDIX B: SOLUTION OF COUPLED INTEGRAL EQUATIONS 



Equation (10) and (11) in the thesis give the two coupled partial 
differential equation that must be solved. Equation (12) gave the form 
that equation (10) and (11) must be in to apply the solution given in 
equation (17). Equation (10) gives the form for the z-directed wire as 



Bz 



4'V r)I + k 2 A z (r) - ^(r) - 



( 10 ) 



Notice that A^(r) i n cylindrical coordinates depends on (R,<j>,z) where 
R is the radius vector in these coordinates. If one fixes R and <J>, then 
A (r) will be a function of z alone along any line parallel to the z-axis . 
Equation (10) then becomes 



Bz 



B rA / N , 2 / v jk P . N B r B , N 

— 7 t[A (z) + k A (z) = L (z) - — [—A (z) ] 

2 z z w z v Bz By y 



( 1 ) 



Equation (12) in the thesis is 



d 2 f(x) 

dx 2 



+ k 2 f(x) = U(x) 



5V(x) 

3x 



( 12 ) 



Examining equation (12) it is seen that equation (1) may be substituted 
directly into equation (12) and equation (17) which is the solution to 
equation (12) may be applied to equation (1). Thus it is seen that for 
equation (10) to use the solution given by equation (17) it must be in 
equation (1) form. This means equation (10) solution is good only along 
a line parallel to z-axis. 
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APPENDIX C: REDUCTION OF A DOUBLE INTEGRAL 



In order to show the solution to equation (47) in general, the sol- 
ution will be merely given and then differentiated to show that it does 
indeed give equation (47) even with arbitrary signs and the addition of 
a constant. The solution is 



(^y + —w)e ^ R cos k(B-e)de = • ? 

R R R[R-(e-f) ] 

where 

R = {a 2 + g 2 U) + [e-f (?) ] 2 } 1/2 . 



-jkR 

{ [e-f ]cos(3-e) 



jRsin k(B-e)} 
( 1 ) 



Now let the solution in general be 



- jkR 



R (R- [e±f (n)±c] ) 
and 



{(e±f±c)cos k(B-e). - jRsin k(g-e)} 



R = {a 2 + g 2 (n) + [a±f (n)±c] 2 } 1//2 , 



( 2 ) 



(3) 



where g and f are functions of the parameter n. 
Now taking the derivative of the function 



d “jkR 

-T-{— n (e±f±c) cos k(B-e) - jR sin k(B-e)]} = 

R -R(e±f±c) 

, "jkR 

[(e±f±c)cos k(B-e) - jR sin k(B-e)] — [— » — =■] + 

R -R(e±f±c) 



[- 



-jkR 

e J d 



R -R(e±f±c) 



j] [(e±f±c)cos k(B-e) - jR sin k(B-e)] 



(4) 
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Working on the first derivative in equation (4) gives 



e 3 2 -ik(> + f+c)e 

6 ] = [ { R -R(e± f ± c) } } - 



dC R 3 -R(e±f±c) 2 



R 



e -jkR | 3R 2 _ ( A ±J _ ± c)_ _ 2R(£±f±c) _ (e±f±c) 2 (e± f±c) j ] 
R R 



3 2 2 

[R -R(e±f±c) ] 



This reduces to give for the first part of equation (4) 



r e ^k^(e±f±c) 



fziii _ i } 

*> R 



R 3 -R(c±f±c) 2 



-] f(e±f±c) cos k(3-c) - jR sin k(g-e)] 



(5) 



Similarly for the second derivative in equation (4) gives the following 
result 



[(e±f±c)k sin k(g-e) + cos k(g-e) + jkR cos K(B-e) 



j sin k(g-e) (e±f±c) n ^ 
R 1 * 



-jkR 



R 3 -R(e±f±c) 2 



( 6 ) 



Adding equations (5) and (6) and factoring out the common elements re- 
sults in 

-jkR 

-a 6 2 t( e± f ±c )k sin k(g-e) + cos k(g-e) 4- jkR cos k($-e) - (7) 

R -R(e±f±c) 



j(e±f±c)sin k(g-e) 
R 



+ 



(e±f±c) (-jk-R)cos k(g-e ) 
R 



(g±f±c) (-jk-R) (-jR sin k[g-e]) , 
R J * 



Now start cancelling common terms 

(e±f±c)k sin k(g-e) - (e±f±c)k sin k(g-e) = 0 
j (e±f ±c)sin k(g-e) j (e±f±c)sin k(g-e) 

_ = C 



(8a) 

(8b) 
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which gives as a result after some algebra 



e -jkR 

~3 2 

R -R(e±f±c) 



^co s k (e_ e )_ { [R 2 _ (e±f±c) 2 Hl + jk R]}) 



R 



(9) 



Cancelling common factor, one has 



e 



-jkR 



R 



R 



-jkR 

e cos 



(3-e) A + 4] 

R R 



( 10 ) 



which is equation (47) in the thesis. 

Thus equation (2) and (3) are required functions in general and 
equation (1) is the necessary function when there is no arbitrary con- 
stant or arbitrary 7 signs in the integrand. 
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